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Abstract 

The origin of both the diffuse high-latitude MeV gamma-ray emission and the 511 keV line 
flux from the Galactic bulge are uncertain. Previous studies have invoked dark matter physics 
to independently explain these observations, though as yet none has been able to explain both 
of these emissions within the well-motivated framework of Weakly-Interacting Massive Particles 
(WIMPs). Here we use an unstable WIMP dark matter model to show that it is in fact possible 
to simultaneously reconcile both of these observations, and in the process show a remarkable 
coincidence: decaying dark matter with MeV mass splittings can explain both observations if 
positrons and photons are produced with similar branching fractions. We illustrate this idea with 
an unstable branon, which is a standard WIMP dark matter candidate appearing in brane world 
models with large extra dimensions. We show that because branons decay via three-body final 
states, they are additionally unconstrained by searches for Galactic MeV gamma-ray lines. As a 
result, such unstable long- lifetime dark matter particles provide novel and distinct signatures that 
can be tested by future observations of MeV gamma-rays. 

PACS numbers: 95.35.+d, ll.10.Kk, 12.60.-i, 98.80.Cq 
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I. INTRODUCTION 



The existence of dark matter is well-established, yet its identity remains elusive. Stan- 
dard dark matter candidates include Weakly- Interacting Massive Particles (WIMPs), which 
have mass ~ 0.1 — 1 TeV, and arise from independent attempts in particle physics to under- 
stand the mechanism of electroweak symmetry breaking. Well-studied signatures of WIMPs 
include elastic scattering off nucleons in underground laboratories, missing energy signals at 
colliders, andparticle production via self-annihilation in Galactic and extragalactic sources 

(e.g. BaaaB). 

These standard signatures, however, do not provide the only means to test the properties 
of dark matter. As an additional example, one may also consider invoking dark matter 
physics to explain anomalous signatures in observed photon or cosmic ray emission spectra, 
even at energies seemingly far removed from those associated with weak-scale physics. In 
fact, in recent years, dark matter models have been constructed to explain astrophysical 
particle production over a variety of energy scales: these include ultra-high energy cosmic 
rays [6], 511 keV line emission from the Galactic bulge 0, @, @, EH HI, H, E3], or the diffuse 
MeV gamma-ray background Il5l |. Though of course all of these anomalies will not be 
due to exotic dark matter physics, the observed emissions can provide strong constraints 
on well-motivated dark matter models, and further they may help identify the interesting 
regions of parameter space for a given dark matter model. 

In this paper we focus on two of the aforementioned anomalies: the diffuse MeV gamma- 
ray background and the 511 keV line flux from the Galactic bulge. As we discuss below, 
known "astrophysical" sources have both spectral shapes and rates that are unable to ac- 
count for these observed emissions. The lack of well- motivated sources, as well as the 
similarities of energy scales, has given rise to speculation that both of these anomalies can 
be explained within the context of a single dark matter model T^, albiet at the cost of 
abandoning the well-motivated WIMP framework. Remaining within the confines of WIMP 
models, these emissions have been separately reconciled by invoking unstable WIMPs with 
nearly de gen erate, ~ MeV scale mass splittings: the 511 keV line flux has been studied in 



Refs. [12], [13|, and the diffuse MeV photon background has been studied in Refs. [17|, [18 



The goal of this paper is to show that, remarkably, both the 511 keV line emission and the 
diffuse MeV gamma-ray spectrum can be explained in the context of a decaying WIMP 
model, characterized by a lifetime of 10 20 s and near equal branching fraction to photons 
and electrons. 

Generically, we focus on a scenario in which the WIMP mass spectrum is highly degen- 
erate, characterized by ~ MeV mass splittings between the next-to-lightest particle (NLP) 
and lightest particle (LP). Both WIMPs freeze out under the standard conditions in the early 
Universe, and one of them is unstable with a lifetime in excess of the Hubble time. The NLP 
decays to three-body final states such as NLP — > LP +7 + 7, NLP — > LP +e + + e~ , NLP 
— > LP +9 + v. Two-body decays are assumed to be either highly suppressed or forbidden. 
Phenomenological consequences of two-body decays with similar lifetimes, and implications 
for the diffuse MeV photon background, were introduced in Ref. (l7| . and further discussed 
in Ref. 19| . Additionally, prospects for detecting long- lifetime decaying dark matter with 
the GeV gamma-ray background were considered in Ref. [201 ] . 

As a specific implementation of the above idea, we consider the brane world scenario 
(BWS), which has become one of the most popular extensions to the Standard Model (SM). 
In the BWS, particles are confined to live on a three-dimensional brane embedded in a higher 
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dimensional (D = 4 + N) space-time, while the gravitational interaction has access to the 
entire bulk space. The fundamental scale of gravity in D dimensions, Md, can be lower 
than the Planck scale, Mp. In the original proposal jU, 22|, the main aim was to address 
the hierarchy problem, and for that reason the value of Mp> was taken to be around the 
electroweak scale. However, brane cosmology models have also been proposed in which Mjj 
is much larger than the TeV scale 23|, |24|. In this paper, we consider a general BWS with 
arbitrary fundamental scale M^; since we neglect gravitational effects, our results do not 
depend on this scale once sufficiently high. 

In general, the existence of extra dimensions is responsible for the appearence of new 
fields on the brane. On one hand, we have the tower of Kaluza-Klein (KK) modes of fields 
propagating in the bulk space, i.e. the gravitons. On the other, since the brane has a 
finite tension, / 4 , its fluctuations will be parametrized by ir a fields called branons. When 
traslational invariance in the bulk space is an exact symmetry, these fields can be understood 
as the massless Goldstone bosons arising from the spontaneous breaking of that symmetry 
induced by the presence of the brane 25, 26[. However, in the most general case, translational 
invariance will be explicitly broken and therefore we expect branons to be massive fields. 
When branons are properly taken into account, the coupling of the SM particles to any bulk 
field is exponentially suppressed by a factor exp[—Mj CK M D /( y 87i 2 f 4 )}, where Mkk is the 
mass of the corresponding KK mode [27, 28]. As a consequence, if / <C Md, the KK modes 
decouple from the SM particles. Therefore, for flexible enough branes, the only relevant 
degrees of freedom at low energies in the BWS are the SM particles and branons. 
The potential signatures of branons at colliders have been considered in Refs. 




3JJ, and astrophysical and cosmological implications have been studied in Refs. 
Moreover in Ref. [33[ the possibility that massive branons could account for the observed 
dark matter of the Universe was studied in detail (see also (HJ). Here, for the first time, 
we consider branon phenomenology with MeV mass splittings and lifetimes in excess of the 
age of the Universe. Due to their universal coupling to the SM, decays of unstable branons 
produce electrons and photons at the same rate in decays, leading to the aforementioned 
consequences for diffuse MeV gamma-rays and the 511 keV line flux. 

This paper is organized as follows. In section [Til we review the observations of the cosmic 
and Galactic gamma-ray backgrounds and the 511 keV signal. Section II II I we summarize 
the salient features of the brane world model, and in section IIVI we present the resulting 
astrophysical signatures. Finally, in sections [V] and ED we discuss other possible signatures 
and recap our main conclusions. 



II. GAMMA-RAY OBSERVATIONS 

In this section, we introduce and discuss the gamma-ray emission spectra we will analyze 
within the framework of decaying dark matter. We focus on two specific observations: the 
511 keV photon line flux from the Galactic center and the high latitude isotropic diffuse MeV 
photon emission. We also, for completeness, discuss the diffuse MeV gamma-ray emission 
from the Galactic center: in the discussion section below we show how these observations 
pertain to our model constraints. 
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A. 511 keV line flux from the Galactic center 



The SPI spectrometer on the INTEGRAL (International Gamma-ray Astrophysics Lab- 
oratory) satellite has measured a 511 keV line emission of 1.05 ± 0.06 x 10~ 3 photons cm -2 



s from the Galactic bulge [34j, confirming earlier measurements [35j. The emission region 
is observed to be approximately spherically-symmetric about the Galactic bulge, with a full 
width half maximum (FWHM) ~ 8°. There is a very low level detection from the Galactic 
disk, ~ 4<t, compared to the 50cx detection from the bulge. The 511 keV line flux is con- 
sistent with an e + e~ annihilation spectrum, fit by the sum of three distinct components: a 
narrow and broad line flux, both centered at 511 keV, and a continuum spectrum extending 
to energies less than 511 keV. The narrow line flux arises both from the direct annihilation 
of thermalized positrons into two photons or through para-positronium formation primarily 
in the cold and warm phase of the inter-stellar medium (ISM). The broad line flux, which 
has a width of 5.4 ± 1.2 keV FWHM, arises from the annihilation of para-positronium in 
flight in the warm and neutral phase of the ISM. The < 511 keV continuum emission arises 
from the annihilations in the ortho-positronium state. The number of 511 keV gamma-rays 



produced per positron is 2(1 — 3p/4), where p is the positronium fraction |36fl. B y detailed 
fitting of the spectrum to annihilation in different phases of the ISM, Ref. |37j] concludes 
that p = O.935i5;0. A model independent fitting of the spectra to a broad line, narrow 
line, continuum spectra, and Galactic diffuse component concludes that p = 0.967 ± 0.022. 
A precise determination of the positronium fraction thus ultimately depends on both the 
nature of the source of the positrons as well as the temperature of the medium in which 
they annihilate. 

The source of the Galactic positrons is uncertain. The fact that the annihilation takes 
place primarily in the warm neutral and warm ionized medium implies that the sources of 
the positrons are diffu sely distributed, and that the initial kinetic energy of the positrons is 



less than a few MeV [361 . |38| . The sources of the positrons are likely contained within the 
observed emission region; the propagation distance from creation to annihilation is at most 
of order ~ 100 pc. Type la supernovae (SNIa) are a candidate for the source of positrons, 
however recent estimates of the escape fraction of positrons from SNIa indicate that they 
cannot account for the entire 511 keV emission 33]. Several other astrophysical sources 
have been proposed jio|, 41, 42|, however none of these sources seem adequate to produce 
the intensity and spatial distribution of the observed line flux. 



B. Isotropic diffuse MeV gamma-ray background at high Galactic latitude 

COMPTEL (the Compton Imaging Telescope) and SMM (the Solar Maximum Mission) 
have measured an isotropic gamma-ray background at high Galactic latitudes (> 10°) over 
the energy ranges 0.8-30 MeV @ and 0.3-7 MeV 0, respectively. More recently, INTE- 



GRAL has measured a diffuse photon background over the energy range 5 — 100 keV [45j. We 



follow the notation of Ref. [19J and refer to these backgrounds as the isotropic diffuse photon 
background (iDPB), as it may include contributions from both Galactic and extragalactic 
sources. At ~ MeV energies, the iDPB analysis is hindered by instrumental and cosmic-ray 
backgrounds, which must be carefully subtracted to reveal the underlying signal. Over the 
COMPTEL and SMM energy ranges, which will be important for our analysis below, the 
observed spectrum is observed to fall like a power law, with dN/dE ~ E~ 2A 

Below energies of a few hundred keV, normal active galactic nuclei (AGN) are able to 
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explain the mean flux of the cosmic X-ray background [46] . A rare population of beamed 
AGN, or blazars, provide an important contribution to the iDPB for energies ^ 10 MeV 47| . 
though there is still room for other sources at these energies |4^ |. The iDPB is observed to 
smoothly transition in between these two energies, however, in the range 1 MeV < E 7 < 
5 MeV, no astrophysical source can account for the observed iDPB. Blazars are observed to 
have a spectral cut-off ~ 10 MeV, and also only a few objects have been detected below this 



energy [49]. SNIa contribute below ~ 5 MeV, but they also cannot account for the entire 
spectrum 50j,|5l|. Recent modeling shows that nonthermal relativistic electrons which alter 
the AGN spectra for energies ^ 1 MeV may account for this excess emission, though a 
detailed understanding of the iDPB at all energies ^ MeV will require matching the mean 



flux and the angular distribution of the sources 52| 



C. Diffuse MeV gamma-rays from the Galactic center 

COMPTEL has determined the flux spectrum of diffuse gamma-rays from the Galactic 



center region over the energy regime 1 — 20 MeV [53j. These fluxes have been averaged 
over a latitude of \l\ < 30° and longitude \b\ < 5°, with high latitudes being used to define 
the zero flux level. It is currently unclear whether the 1 — 20 MeV spectrum is a result of 
diffuse or point source emission. For energies > 100 MeV, the gamma-ray spectrum is likely 
produced by both nucleon interactions with interstellar gas via neutral pion production and 
electrons via inverse compton scattering. For energies < 100 keV, point sources dominate 



the gamma-ray spectrum 54]. However, an inverse compton spectrum that matches the 



diffuse gamma-ray spectrum at ~ 100 MeV can account for at most 50% of the emission 



between 1-20 MeV 55 



In addition to the diffuse measurement from COMPTEL, INTEGRAL has performed a 
search for gamma-ray lines originating within 13° from the Galactic center. Two lines were 
recovered: the aforementioned line at 511 keV and an additional line at 1809 keV. The 1809 
keV line originates from the hydrostatic nucleosynthesis of radioactive elements in the cores 
of massive stars, where long-lived isotopes such as 26 A1 are carried to the ISM, after which 
they subsequently decay mono-energetically. The 1809 keV line flux has been detected by 
INTEGRAL at a level of ~ 10~ 4 cm -2 s _1 . Other than these two gamma-ray lines, no other 
lines were detected up to upper flux limits of 10 _5 -10~ 2 cm -2 s -1 , depending on line width, 
energy, and exposure. 

III. BRANON PHENOMENOLOGY 

Having introduced and discussed the gamma-ray emission spectra we consider, in this 
section we will briefly review the main properties of massive brane fluctuations (see Refs. [2H 
for a more detailed description). Focusing in particular on models that produce the 



observed abundance of dark matter, we present the formulae for determining the decay 
widths of long-lifetime branons to photons, electrons, and neutrinos. 

A. Br anon overview 

We consider a typical brane model in large extra dimensions. The four-dimensional space- 
time, M 4 , is embedded in a D-dimensional bulk space which, for simplicity, it is assumed to 
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factorize as Mr, = M4 x B. The extra space B is a given N- dimensional compact manifold, 
so that D = 4 + N. The brane lies along M4 and we neglect its contribution to the bulk 
gravitational field. The bulk space Md is endowed with a metric tensor Gmn, which we will 
assume for simplicity is given by 

° MN ~ [ -g' mn {v) ) ■ (1) 

The position of the brane in the bulk can be parametrized as Y M = (x fl ,Y rn (x)), with 
M — 0, ... ,3 + N. We have chosen the bulk coordinates so that the first four are identified 
with the space-time brane coordinates x^. We assume the brane to be created at a certain 
point in B, Y m (x) = Y™, which corresponds to its ground state. We will also assume that 
B is a homogeneous space, so that brane fluctuations can be written in terms of properly 
normalized coordinates in the extra space: n a (x) = f 2 Y a (x), a = l,...,N. The induced 
metric on the brane in its ground state is simply given by the four- dimensional components 
of the bulk space metric, i.e. g^ u = g^ u = G M „. However, when brane excitations are present, 
the induced metric is given by 

9^ = d^d^GMNix, Y(x)) = g^x, Y(x)) - d,Y m d u Y n g' mn (Y(x)) . (2) 

The contribution of branons to the induced metric is then obtained by expanding Equation [2] 



around the ground state [26|, |30|, [56 



g»v = g»v - jiScfld^dyif + jpg^M^^ir + ... (3) 

Branons are the mass eigenstates of the brane fluctuations in the extra-space directions. 
The branon mass matrix M a p is determined by the metric properties of the bulk space and, 
in the absence of a general model for the bulk dynamics, we will consider its elements as 
free parameters (for an explicit construction see Refs. (H3|). Therefore, branons are massless 



only in highly symmetric cases |26|, [30|, [56|, [58 



Since branon fields survive in the limit in which gravity decouples, Md — * 00, branon 



effects can be studied independent of gravity 59] . We will work in the thin brane limit and 



assume that the brane dynamics can be described by a low-energy effective action derived 
from the Nambu-Goto action j26|. Also, branon couplings to the SM fields can be obtained 
from the SM action defined on a curved background given by the induced metric Equation^ 
and expanding in branon fields. Thus the complete action, up to second order in ir fields, 
contains the SM terms, the kinetic term for the branons and the interaction terms between 



the SM particles and the branons [25j, |26|, |30|, [56 



It is interesting to note that under a parity transformation on the brane, the branon field 
changes sign if the number of spatial dimensions of the brane is odd, whereas it remains 
unchanged for even dimensions. Accordingly, branons on a 3-brane are pseudoscalar parti- 
cles. This fact, in addition to the geometrical origin of the action, implies that terms in the 
effective Lagrangian with an odd number of branons are forbidden. It means that branons 
are stable. Strictly speaking, this argument only ensures the stability of the lightest branon 
(7T 1 ), the rest of branons can decay to the lightest one, producing SM particles. These decays 
are very suppressed if the branon spectrum is very degenerate. 

In order to analyze such properties we will study the simplest model containing an unsta- 
ble branon. We consider a model where we have two branons, implying that the number of 
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extra dimensions has to be two or greater: the lightest one is defined as tt 1 , and the heaviest 
one as 7f h . We will assume very similar masses for them: AM = M%h — M%\ -C M ~ M%h ~ 
Mfti. Using these energ y eig enstates, the SM-branon low-energy effective Lagrangian may 
be written as 25|, |26|, [30j, |5£ 



*— "K- 



(4) 



Here a((3) = h, I, and T£ v is the standard energy-momentum tensor of the particle evalu- 
ated in the background metric: 



T, 



ere, 



+ 2^ 



(5) 



M, whereas the interaction 



The mass matrix is diagonal with eigenvalues J4h ~ AL. 
matrix, I a/3 , is exactly the identity at first order 25, 26,[3(J |56j. However, radiative correc- 



tions and higher order terms have a different effect on the interaction eigenstates relative 
to the energy eigenstates, resulting in suppressed but non-zero cross interaction parameters, 



l aj3 



1 A /2 
A^/2 1 



(6) 



This implies that, typically, A^ ~ 0.01. For simplicity, we will take this parameter to be 
independent of the particular SM particle, i.e. A^ ~ A for any <p. The main conclusions of 
this paper do not depend on this assumption, although it does introduce an uncertainty of 
order one for a particular A^,. However, for the most part of branon phenomenology, both 
this parameter and AM are negligible. General constraints on the branon parameter space 
we consider are shown in Figured) the constraints include those from HERA, Tevatron, and 
LEP-II (additional bounds from astrophysics and cosmology can be found in 60j). 

Thus A is fundamental to analyze the stability of the heaviest branon. Indeed, if A is 
exactly zero, 7f h is completely stable. If A is non zero, the heaviest branon decays to the 
lightest branon and a Standard Model particle anti-particle pair. For the analysis below, we 
are interested in AM < 10 MeV, which means that the only available decay channels are 
T e + : 7T h — > fc l e + e~ , r 7 : n h — > ff 77, and T u : 7r h — ► n uu. 



B. Thermal abundances 

Branons thus interact in pairs with the SM particles, and the lightest is necessarily 
stable. In addition, their couplings are suppressed by the brane tension f A , which means 
that they could be in general weakly interacting and massive. As a consequence their freeze- 
out temperature can be relatively high, which implies that their relic abundance can be 
cosmo logically significant. 

In order to calculate the thermal relic branon abundance, the standard techniques have 



been used in the case of non-relativistic branons at decoupling [33|, l60|. The evolution of 
the number density n a of branons interacting with SM particles in an expanding Universe 
is given by the Boltzmann equation: 

''" " -3Hn a -(a A v)(n 2 a -(n^f) (7) 



dt 



7 



2500 



2000 



1000 



500 



ACCESSIBLE 



CLIC ^^^^^p 


ACCESSIBLE y ^\^r 
ILC V 




/ EXCLUDED 
/ LEP + TEVATRON 


ACCESSIBLE 
LHC 





500 



1000 

M[GeV] 



1500 



2000 



FIG. 1: The relic branon abundance for N = 2 in the f — M plane (see Refs. 33|, |6(| for details). 



The lower area is excluded by single-photon processes at LEP-II 30|, [62( together with monojet 
signals at Tevatron-I 0. Prospects for real branon production at future colliders are also shown 



(See Refs. 30j, [63|). The region between the two horizontal dotted lines are preferred by the muon 



anomalous magnetic moment and electroweak precison observables 28|, |61 ] . 



where 

<7 A = X>(7r<V - X) (8) 

X 

is the total annihilation cross section of branons into SM particles X summed over final 
states. The —3Hn a term, with H the Hubble parameter, takes into account the dilution of 
the number density due to the expansion of the Universe. This is the dominant mechanism 
that sets the total number density of branons. In the scenario we are interested in, the 
heaviest branon will also decay to the lightest branon, however this effect is insignificant in 
changing their respective number densities because we focus on lifetimes longer than the age 
of the Universe. Therefore, the relic density is determined only by the parameters M and 
/. Each branon species evolves independently and has exactly the same abundance before 
the decays are effective. 

The relic abundance results for two branons of mass M are shown in Figure [TJ As 
expected, for / and M scales of ~ 100 GeV — 1 TeV, we have the correct amount of total 
non-baryonic dark matter abundance. As shown, these scales are not only in the natural 
region of parameter space, they are also in the favored region for Brookhaven determinations 



of the muon anomalous magnetic moment [28|, [61 
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C. Branon decays 



As we have seen, decays of very degenerate branons can be described at low energies by 
an effective action that depends on four parameters: the branon mass M, the mass splitting 
AM, the cross interaction parameter A, and the brane tension scale /, which suppresses the 
coupling of these new particles with the SM. We now determine the decay widths to the 
relevant SM particles. 

Since branons couple directly to the energy momentum tensor, they do not couple directly 
to a single photon, and thus decays into two-body final states that include a single photon 
are forbidden. Indeed, this is a general property of scalar or pseudo-scalar particles; the 
decay of a spin-zero particle into another spin-zero particle and a photon is forbidden by 
angular momentum conservation. On similar grounds, the decay of a spin-zero particle into 
another spin-zero particle and a fermion is forbidden. Therefore, unless another light spin- 
zero particle is added to the Standard Model, the decays of degenerate scalar multiplets will 
proceed predominantly into three-body final states. 

We focus specifically on the limit M » AM. In this limit, the decay widths can be 
calculated from Equation HI the branon Feynman rules given in Ref. (30|, and substituting 
S a/ 3 for I a p. For photons, we get a differential decay width of 



rfr 7 (e 7 ) _ A 2 4(AM-e 7 ) 3 M 2 



de 1 24 / 8 7t 3 



which implies a total photon decay width of 

2 



1.67 x 10 20 s. 



10 



-2 



2A 



J_ 
M 



4MeV 
AM 



/ 



ITeV 



(9) 



(10) 



If AM > 2m e , the heaviest branon can also decay to the lightest one and an electron-positron 
pair. The spectrum of the outgoing positron in the same limit of M >> AM is: 



dT e +(e e +,m e 
ds P + 



A M 2 



9 ( |7r) - \l ( £ l+ ~ m 1)((AM - e e +) 2 - ml) 

{4(AM -2e e+ ) 2 [(AM - e e+ )e e+ - 16 m 2 ] 
+ m 2 [225 (AM - e e+ )e e+ - 177m 2 ]} . 



(11) 



The kinematic limits for the total energy of the outgoing positron are m e > e e + > AM — m e . 

Finally, we note that in the kinematic limits we study it is also possible to produce 
neutrinos: 



de Vi 



lrfr e+ (e„.,0) 
2 de Vi 



(12) 



where i labels the species of neutrino (we have neglected the neutrino mass). Since diffuse 
neutrino bounds are much weaker than diffuse photon bounds at MeV energies, there are 
likely minimal phenomenological implications of neutrino production in this model (see 
Refs. (641 . 65 ] for recent determinations of neutrino constraints of annihilating and decaying 



dark matter models). The decay widths for all of these channels are shown in Figure [2j 
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FIG. 2: Left: The differential decay widths for 7r h 



e , fi" 77, and tt vis, assuming AM 



3.5 MeV. Right: The total branching fraction for each function of AAf. In both figures, 

for isD, we have accounted for the decays to all three species of neutrinos. 

IV. GALACTIC AND EXTRAGALACTIC SIGNALS 

We now determine the extragalactic and Galactic gamma-ray signals from the decays 
of unstable 7r h, s. Although the resulting fluxes of course depend on the specific model 
parameters and couplings, the formalism we present here can be used for any similar model 
that decays via three-body final states. In what follows, the mass of the NLP is defined as 
M, the mass splitting between the NLP and LP is AM, and the lifetime of the NLP into a 
given branching ratio, i, as 1/r* = T t = J dEdTjdE. 

A. Galactic 511 keV flux 

We begin by considering the 511 keV line emission resulting from positron injection in the 
Galactic bulge. Given the small propagation distance of the positrons as described above, 
we assume that the location of the e + e~ annihilation and 511 keV photon creation faithfully 
tracks the dark matter halo distribution. Defining s as the distance from the Sun to any 
point in the halo, \& as the angle between the Galactic center and any point in the halo, and 
D = 8.5 kpc as the distance from the Sun to the Galactic center, the 511 keV line intensity 
from decays is 



where r 2 (s, = D 2 + s 2 — 2Ds cos\l/ , and the solid angle is defined as AQ = 2n(l — 
cos^). The number density of dark matter particles is p(r)/M, where p(r) is the density of 
dark matter in the halo of the Milky Way, which we assume to be spherically symmetric. 
Consequently, the total number of 511 keV photons produced per unit time can be estimated 
as dn 5 u/dt = 2(1 — 3p/4)n^hT e + = (1 — 3p/4) pT e +/M, where we are assuming that n h 
accounts for half of the dark matter (the other half is in form of if 1 ), that the positrons 
are stopped before annihilation, and that this annihilation takes place through positronium 




(13) 
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formation p fraction of the time. As we have discussed above, p ~ 0.94. 

Both the total flux within the field-of-view and the angular distribution of the flux depend 
crucially on the shape of the dark matter density profile. Numerical simulations with only 
dark matter have shown that the halos consisting of cold dark matter particles have a 
Navarro-Frenk- White (NFW) type density profile, p(r) = p /(r/r )/(l + r / r o) 2 ? where po is 
the scale density and r is the scale radius [66]. The NFW profile, however, does not account 
for the effect of baryonic physics, which is expected to be important in setting the central 
density of Milky Way-type galaxies, due to the relatively large baryonic contribution in the 
central regions. A full understanding of the halo density profile must adequately account 
for the energy exchange processes between the baryons and dark matter, which is currently 
lacking. These interactions are particularly important in the interior regions, and it is 
currently not understood if these interactions steepen the central density of the dark matter 
halo or flatten it out. From an observational perspective, there is also a wide uncertainty 
in the density profile and mass model of the Milky Way halo. For example, Ref. 67| argues 
that the large number of microlensing constraints towards the Galactic bulge is inconsistent 
with central slopes steeper than r~ 0A . Ref. [68| apply a variety of observational constraints 
to an adiabatically- contracted NFW model, and find that NFW-like central slopes provide 
an excellent fit to the data. 

To allow for the greatest flexibility, we model the density profile of the Milky Way in the 
form of 

P(r) = nr^T' ( 14 ) 

(r/r o y[l + (r/r o rf^ )/a 

where po is the scale radius, r is the scale density. The combination of 7 and (3 set the 
inner and outer slopes of the halo profile, respectively, while a controls the sharpness of 
the transition between the inner and outer slopes. We determine the shape parameters 
a, (3, and 7, as well as the normalization parameters po and r by fitting to both the 
angular distribution of the INTEGRAL signal and the constraints on the Milky Way halo. 
The parameters po and 7 are primarily determined by the INTEGRAL data, whereas the 
remaining parameters are fixed by the total mass and local density of the dark matter halo. 
The intensity of the INTEGRAL signal is so low at high longitudes that it is difficult to 
discriminate from the instrumental background, and thus fix the parameters that do not 
directly effect the central region of the halo. 

In Figure [3], we show an example halo profile that fits the angular distribution of the 511 
keV intensity, and is normalized to match the total flux of 1.05 x 10~ 3 cm" 2 s _1 within the 
Galactic bulge. For simplicity, we assume that the 511 keV emission is entirely a result of 
dark matter decays, and does not include any contribution from sources discussed above, 
such as SNIa. Here we use AM = 3.5 MeV, = 700 GeV, A = 0.0085, and a brane 
tension of / = 1 TeV; these model parameters naturally give the correct relic dark matter 
abundance, as seen in Figure [TJ The implied parameters describing the halo model are 
Po = 0.12 GeV cm -3 , r = 10 kpc, 7 = 1.5, (3 = 3, a = 8. With these halo model 
parameters, the total dark matter halo mass is in accord, to within statistical uncertainties, 



with the mass determinations from the kinematics of satellite galaxies [69j, [70J • In addition 
to the true signal, which is shown as a solid line in Figure [31 we show two curves with angular 
resolutions of 1° and 3°; these are similar to the ~ 3° FWHM angular resolution of SPI. 

In Figure [31 we have focused on matching the normalization and angular distribution 
of the flux within the Galactic center region, where the 511 keV flux is well-determined. 
As an additional check on the viability of the model, we must be sure not to overproduce 
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FIG. 3: The 511 keV intensity from the decays of unstable fr h 's as a function of Galactocentric 
angle. The intensity has been normalized to the total flux of 1.05 x 10~ 3 cm -2 s _1 within the 
Galactic bulge, and the dark matter halo is described by Equation O with a = 8, (3 = 3, and 
7 = 1.5. The solid (black) line is the true signal, the dashed (red) line is the signal smoothed with 
a 1° angular resolution, and the dotted (blue) curve is for a 3° angular resolution. We have used 
a positronium fraction p = 0.94, and the branon model parameters AM = 3.5 MeV, / = 1 TeV, 
and Mjfh = 700 GeV. The angular distribution of the INTEGRAL signal is spherically- symmetric 
about the Galactic bulge, with a FWHM of 8°. 

the observed 511 keV emission from high longitudes in the Galactic disk. We can examine 
constraints on a potential high longitude disk flux by separately considering both emission 
from the Galactic center region and large-scale extended diffuse emission. Regarding the 
Galactic center region emission, we use the most recent INTEGRAL results, which shows 
that the 511 keV bulge-to-disk flux ratio is ~ 1—3 and that the intensity is reduced by a factor 
~ 6 at latitudes and longitudes of 10° relative to the Galactic center [34J. As seen in Figure [3l 
the chosen profile satisfies these constraints. Regarding the large scale (> 20°) diffuse 511 
keV emission, we consider the results of Ref. [711 ] . Using the INTEGRAL effective area of 
75 cm 2 at 511 keV, and smoothing the solid curve in Figure [3] with a gaussian comparable 



to the ~ 20° FWHM spatial scales resolved in Ref. [7l|, we find a 511 keV rate of ~ 0.06 s 1 



at / = 0°, and ~ 0.01 s" 1 at / = 30°. These rates were determined by subtracting the flux 



at b = 20°, and are consistent with the rates determined in Ref. [71j. Thus we find that the 
halo profile used in Figure [3] is consistent with 511 keV emission constraints at high Galactic 
longitude. 

Before moving on to consider the direct production of gamma-rays in 7r h decays, we make 
one additional comment regarding the assumed value of the positronium fraction. Motivated 
by the detailed fits by INTEGRAL, for all of the above results we have taken the positronium 
fraction to be p = 0.94. However, the precise determination of p is model dependent, as 
discussed above. If the decays we consider do indeed account for a large fraction of the 511 
keV signal, a redetermination of p will be required by explicitly fitting both the Galactic and 
extragalactic fluxes. Inclusion of the latter flux will be important, even in the direction of 
the Galactic center, given that the extragalactic and Galactic fluxes are similar to within an 
order of magnitude, depending specifically on the shape of the halo model. The similarity 
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of extragalactic and Galactic fluxes is unique to the models we consider; in the case of 
dark matter annihilation, for example, the extragalactic flux is typically several orders of 
magnitude below the flux from any direction in the Galaxy itself. 



B. High latitude gamma-ray flux 



We now turn to determining the contribution to the high latitude iDPB from three-body 
decays to gamma-rays, both from the Galactic halo and from extragalactic sources. First 
considering extragalactic sources, these decays will be visible only if they occur in the late 
Universe, in the matter or vacuum dominated eras. Assuming a smooth distribution of 
dark matter in the Universe (modifications to this assumption will be discussed below), the 
differential gamma-ray flux from a general three-body decay is 



rf$ _ n 7 r*> N(t)dT 
dE~ An Jo aV de 



7 



(15) 



7 



Here n 7 is the number of photons produced in a single decay; a is the scale factor of the 
Universe; to ~ 4.3 x 10 17 s is the age of the Universe; N(t) = N m e~ t / T , where N m is the 
number of decaying particles at freeze-out; and Vq is the present volume of the Universe. The 
relation between the produced energy, e 7 , and observed energy, E 7 , is given by de 1 /dE 1 = 
(1 + z) = a -1 . In Equation [TSJ we have assumed that gamma-rays do not get attenuated on 
cosmological gamma-ray backgrounds from their redshift of production. We find this to be 
an excellent approximation over the energies and redshifts we are considering [72|]. In the 
BWS, n 7 = 2, however this factor gets cancelled because 1/2 of the dark matter consists of 
unstable 7r h 's. 

We are interested in decays such that r ^> £ . In this case, we can write t ~ P(a), where 



P(a) 



2 (in V^a^ 3 + V^m + ^ao 3 — In v/OjVf 



The differential photon flux coming from the decays may be written as 



n 7 iV m f 1 



dE„ 



e -P(a)/r dT 

da 



An V JE-y/AM aQ(a) de 



where we have defined 



(16) 



(17) 



da 
~dt 



aH ~ H n 



in 



M 



+ Q A a 2 = Q(a). 



(18) 



For lifetimes much longer than the age of the Universe, the exponential factor in Equation [T71 
is negligible. We neglect the radiation content, ~ 0, as well as the curvature term, k ~ 0. 
Though the extragalactic flux is dependent on cosmological parameters [13], here we choose 
to fix the parameters to match the observed concordance cosmological model [73|. 

The high latitude iDPB also receives a significant contribution from decays in the halo of 
the Galaxy. The flux from Galactic halo dark matter decays directly into photons can also be 
calculated using Equation [131 accounting for the fact that the flux comes from high latitudes, 
^ 10°, and averaging the flux over a ~ 1CT 3 sr field of view. This latter quantity corresponds 
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FIG. 4: The high latitude isotropic diffuse photon background in the MeV regime, as measured by 



COMPTEL (circles) 43]. The short-dashed curve is the contribution from extragalactic decays, 
the long-dashed curve is the contribution from decays in the Galactic halo, and the solid curve is 
the sum of both extragalactic and Galactic contributions. The branon model parameters are the 
same as those in Figure Two additional standard contributions to the iDPB are shown: Type la 
Supernovae (SNIa; dotted) and blazars (dot-dashed). The blazar background has been normalized 



to the gamma-ray background at energies ^ 10 MeV 49(], while the normalization of the SNIa 
spectrum is taken from Ref. [icf . 

to the field of view characteristic of MeV gamma-ray instruments such as COMPTEL. With 
BWS parameters fixed by the 511 keV signal, in Figure H] we show the contribution to the 
iDPB from high latitude Galactic and extragalactic 7r h decays. Again, unstable 7r h, s account 
for 1/2 of the relic abundance of dark matter. In the calculation of the observed spectrum, 
we have not accounted for finite detector energy resolution; including these effects will likely 
further smooth out the endpoints of the spectrum, and bring the shape of the spectrum into 
even better agreement with the approximate power-law spectrum observed by COMPTEL. 



V. DISCUSSION 

Above we have shown that, assuming both an abundance of non-baryonic dark matter 
that matches cosmological observations and WIMP scale masses, within the BWS two pa- 
rameters are required to independently match the 511 keV line flux and the ~ MeV iDPB. 
We now discuss independent constraints on this model as well as potential implications of 
these results. 



A. Gamma-rays from the Galactic center 

In section IIVI we have focused on the high latitude Galactic gamma-ray background. 
There will additionally of course be a significant contribution to the low-latitude gamma- 
ray background from n h decays, in particular in the direction the Galactic center. We must 
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be sure that the model parameters we consider above do not violate the COMPTEL flux 
measurements in this direction. 

We can make a simple estimate of the gamma-ray flux in the direction of the Galactic 
center by considering the branching ratios determined in section IIHI Using these branching 
ratios, in combination with the number of gamma-rays produced per e + e~ annihilation 
(section [II]), the ratio of the total number of photons produced at 511 keV to the total 
number of photons produced directly in the two photon decays for our model parameters is 
simply 

^ (l - -p) ~ 0.9. (19) 



r 7 v a 1 



In general, we can expect a ratio of order one for any model that has similar branching 
ratios to photons and positrons. Thus we expect a continuous gamma-ray flux < 10~ 2 
cm -2 s _1 sr _1 at ~ MeV energies from the Galactic center. Interestingly, COMPTEL has 
measured a diffuse flux from the Galactic center at this order of magnitude in the energy 



range 1 — 3 MeV [55] . The COMPTEL flux is seen to vary slowly with energy, and is reduced 
by a factor ~ 2 at a Galactic latitude of ~ 10°. We find that, for the dark matter halo profile 
considered above, the diffuse MeV gamma-ray emission from 7? h decays is consistent with 
the COMPTEL flux both at the Galactic center and at higher latitudes. Since, as discussed 
above, the source of MeV gamma-rays from the Galactic center is uncertain, it is interesting 
to consider dark matter decays as the source of this emission. 

In addition to the diffuse Galactic gamma-rays, INTEGRAL has placed strong constraints 
on the presence of anomalous gamma-ray lines in the energy regimes ~ 10 — 8000 keV within 



13° from the Galactic center 711 ] . For example, at energies ~ MeV, lines of width ~ keV 
are constrained to have a flux < 10~ 4 cm -2 s _1 , and "lines" of width ~ MeV are constrained 
to have a flux % 10~ 3 cm~ 2 s _1 . Given the width of the gamma-ray spectra we consider (see 
Figure [2]), the latter flux limits are relevant for 7f h decays. As a result, the INTEGRAL line 
constraints are no more stringent than the COMPTEL constraints on diffuse MeV gamma- 
ray emission. 



B. Assumption of smooth halo and angular distribution 

The extragalactic flux from dark matter decays depends on the distribution of dark matter 
halos in the Universe, while the Galactic flux depends on the distribution of dark matter 
within the Milky Way halo. In principle, both of these fluxes also depend on the presence of 



'substructure" within dark matter halos 74 . While we have assumed that the distribution 



of dark matter is smooth both in the Universe and in the Galactic halo, we can gain some 
insight as to the effects of clumping by considering the simplifying limit in which all of the 
substructures are at a common mass, the internal density of the substructures is constant, 
and the radial distribution of substructures in the halo follows the same density profile as 
the smooth dark matter halo. In reality, it is likely that the two radial distributions differ, 
which may affect our simple estimate. Under these assumptions, we find that, in a given 
direction, the increase in the flux over the smooth halo is simply 1 + /, where / is the fraction 
of the halo contained in substructures. Typically, / ~ 0.01 74J. This small increase simply 
quantifies the probability that any one of these substructures is within a given field of view. 

Regarding the determination of the extragalactic gamma-ray flux from 7f h decays, we have 
assumed a smooth distribution of 7f h 's in the Universe, and have neglected the clumping in 
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dark matter halos. Given the above considerations, the clumping in halos will likely have 
minimal effect on the mean flux signal, though the situation will be different when consider- 
ing the angular distribution of gamma-rays. The angular distribution depends not only on 
the parameters describing the branon model, but also on the cosmological evolution of dark 
matter halos as a function of redshift. Similar to the case of standard annihilating neutralino 
WIMP dark matter, the angular distribution of gamma-rays from decays is expected to be 
distinct from astrophysical sources such as AGN and super novae, and may provide a way 



to discriminate between the sources of the iDPB [75j, [76J. We consider this analysis in a 
forthcoming study. 



VI. CONCLUSIONS 

In this paper, we have analyzed the general possibility that decaying dark matter may be 
the source of both the 511 keV gamma-ray line from the Galactic center and the diffuse high 
latitude photon background at energies between 1 and 5 MeV. We have introduced a new 
scenario where a dark matter particle is nearly degenerate with a lighter daughter, with mass 
splitting ~ MeV and lifetime ~ 10 3 times the present age of the Universe. The decays of the 
dark matter are into three body final states: a Standard Model particle anti-particle pair 
and a daughter dark matter particle. Remarkably, we find that the above mass splitting 
and lifetime are exactly that which is required to explain both of the above gamma-ray 
emissions, provided there is similar branching ratios to electrons and photons. 

We have illustrated this idea with a concrete model of brane-world dark matter, and 
we have shown that a standard WIMP motivated within the brane-world scenario, the 
branon, can explain both observations with natural values required to obtain the correct 
dark matter relic abundance. Using the effective low-energy Lagrangian for a model with 
two branons interacting with the Standard Model fields, we have computed the decay rates of 
the heaviest branon into the lighter one and pairs of Standard Model particles, in particular 
focusing on electron-positron pairs and photons. We find that the branon model not only 
produces positrons and gamma-rays at the required rates, it also provides unique and novel 
phenomenology that will be tested by future MeV gamma-ray observatories. 
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